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Abstract. 

We study the position of a biased tracer particle (TP) in a bath of hardcore particles 
moving on a lattice of arbitrary dimension and in contact with a reservoir. Starting 
from the master equation satisfied by the joint probability of the TP position and the 
bath configuration and resorting to a mean-field-type approximation, we presented a 
computation of the fluctuations of the TP position in a previous publication [Phys. 
Rev. E 87, 032164 (2013)]. Counterintuitively, on a one-dimensional lattice, the 
diffusion coefficient of the TP was shown to be a non-monotonic function of the 
density of bath particles, and reaches a maximum for a nonzero value of the density. 
Here, we: (i) give the details of this computation and ofi'er a physical insight into the 
understanding of the non-monotonicity of the diffusion coefficient; (ii) extend the mean- 
field-type approximation to decouple higher-order correlation functions, and obtain the 
evolution equation satisfied by the cumulant generating function of the position of the 
TP, valid in any space dimension. In the particular case of a one-dimensional lattice, 
we solve this equation and obtain the probability distribution of the TP position. 
We show that the position rescaled by its fluctuations is asymptotically distributed 
accordingly to a Gaussian distribution in the long-time limit. 
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1. Introduction 

1.1. Context 

Studying the dynamics of an active particle or a particle submitted to an external force 
travelling in a crowded environment is a frequent problem in physics and in biology. 
Different examples are found in biophysics, when one considers molecular motors, motile 
living cells or bacteria BEl. or in the study of biased intruders in granular systems 
[3] or colloidal suspensions |1]. The determination of the dynamics of this tracer 
particle (TP) is consequently an important question, with different applications. In 
particular, new experimental tools allow the study of a medium using a microscopic 
probe particle submitted to an external force. This field of research, commonly named 
active microrheology, has been applied to different systems in the past decades, among 
which biological cells m, complex fiuids and colloidal suspensions m 

From a theoretical point of view, the difficulty lies in the modeling of the 
environment of the tracer particle (TP), which is constituted of a large number of 
interacting bath particles. In most analytical approaches, the evolution of the position 
of the TP is studied with some effective description of the bath of particles HDl, which 
do not take into account the correlations between the position of the TP and the density 
profiles of the surrounding bath. In the situation where the TP and the bath particles 
have comparable sizes, the response of the probe to the external forcing, and in particular 
its fluctuations, cannot be accounted for correctly with an effective treatment. 

Here, we study the diffusion of a driven tracer in a lattice gas of hardcore parti¬ 
cles. This minimal model explicitly takes into account the bath dynamics: the driven 
TP performs a biased nearest-neighbor random walk, in a bath of particles performing 
symmetric nearest-neighbor random walks, with the restriction that each site is occu¬ 
pied by at most one particle. We assume that the lattice is in contact with a reservoir 
of particles, so that the bath particles present on the lattice may desorb back to the 
reservoir, and particles from the reservoir may adsorb onto vacant lattice sites. This 
so-called Langmuir kinetics is relevant to describe situations where a gas or a vapor is 
brought in contact with a solid surface, on which the gas particles may form an adsorbed 
layer. The transport properties of the adsorbed particles have been shown to control 
many different processes, such as spreading of molecular films on solid surfaces m or 
dewetting [121 US]- The particular case where the Langmuir kinetics is coupled to a 
Totally Asymmetric Exclusion Process was investigated theoretically [HE], and has 
been show to be relevant to describe the directional motion of molecular motors on a 
cytoskeletal filament, with random attachment and detachment of the motors [isiin]. 
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Studying the transport properties of a biased TP in a hardcore lattice gas is actnally 
a complex iV-body problem. In the sitnation where the density of bath particles p is 
very high and where the number of particles on the lattice is conserved, the problem 
can be treated exactly to obtain resnlts at leading order in (1 — p) [TSl UHl 1201 [H]. For 
an arbitrary density of particles, exact results were established concerning the mean 
position of the TP in the one-dimensional situation |22], and concerning the validity of 
the Einstein relation |23]. The sitnation where the lattice is populated by an arbitrary 
density of particles and where it is in contact with a reservoir of particles was addressed 
by resorting to a mean-field-type approximation consisting of the deconpling of relevant 
correlation fnnctions, allowing the computation of the mean position of the TP and of the 
bath density prohles in the long-time limit, in the case of a one-dimensional lattice |21] 
and for lattices of higher dimension |25l EHl EZl EH]- Nnmerical simnlations sampling 
exactly the master eqnation of the problem revealed the accuracy of the decoupling 
approximation in a wide range of parameters. The situation where the bath particles 
are hxed and can appear and disappear with prescribed rates, known as dynamical 
percolation IHH], was also considered [HU] . 

More recently, by extending the decoupling approximation initially proposed to 
stndy the mean position of the TP, its fluctnations have been stndied [HT] . An evolntion 
equation for the flnctnations of the TP on a lattice of arbitrary dimension was obtained. 
This equation was solved explicitly in the case of a one-dimensional lattice, and in the 
stationary limit. The analysis of the solutions revealed a striking feature of the diffusion 
coefficient: in a wide range of parameters, it is shown to be a nonmonotonic fnnction of 
the density of particles on the lattice. Counterintnitively, the diffusion coefficient of the 
TP can then be enhanced by the presence of bath particles in its environment. 

Recently, this nonmonotonicity was also observed in the sitnation where the TP is 
dragged in a bath of soft particles |32l [33], and then appears to be a generic featnre of 
biased intruders in crowded environments. 

1.2. Main results of this paper and overview 

In this paper, we hrst give a detailed computation of the resnlts presented in [31]: we 
establish the evolntion equation of the fluctnations of the TP position in arbitrary 
dimension nnder the decoupling approximation, and solve it in the case of a one¬ 
dimensional lattice. We show that the non-monotonicity of the diffnsion coefficient with 
respect to the density of particles on the lattice is correlated to a non-monotonicity of 
some cross-correlation fnnctions with respect to the distance to the TP. 

The main result of this paper is the following: we go one step further and generalize 
the mean-field-type approximation in order to calculate the cumulant generating 
fnnction of the position of the TP, and therefore its complete probability distribntion. 
Denoting by Xt the position of the TP, by Xt = XfCi its projection along the direction 
of the bias, and defining the cumnlant generating fnnction \k(M; t) = ln(e'“^‘), we obtain 
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in the long-time limit 


'^{u;t) ~ ^{u)t, 

t^OO 


with 


$(m) = ^ (e™" - 1) 

r 


1 - 


(e‘^^*?7x,+ei) 

l^QiuXt ^ 


r ^ 


-1) 


1 - 


(e‘“^*77Xt-ei) 

{QiuXt) 


( 1 ) 


, ( 2 ) 


where p^ is the jump probability of the TP in the direction v of the lattice, r is its 
characteristic jump time, a is the lattice spacing and rjr is the occupation number at site 
r. The evolution equations for the correlation functions Wr{u; t) = (e™^*pxt+r)/(e^“^‘) 
are obtained using an extension of the mean-field-type approximation proposed to study 


the mean and fluctuations of the position of the TP. This result is given by Eqs. (74) 


and (75), where the quantities kr{t) = {rjxt+r) are the density profiles in the reference 


frame of the TP and are the solutions of Eqs. (15) and (16). The equations satisfied by 


the correlation functions Wr{u]t) are solved in the particular case of a one-dimensional 
lattice, and we compute the probability distribution of the position of the TP. We show 
that all the cumulants of the TP position scale as t. 

We also consider the random variable obtained by rescaling the position of the TP 
by \/i. For this random variable, all the cumulants of order greater than 2 vanish in 
the long-time limit, which shows that the position of the TP is asymptotically Gaussian. 


The article is organized as follows. In section 2, we present the model and give 
the master equation governing the joint probability of the position of the TP and of 
the bath particles configuration. In section 3, we give the evolution equations of the 
first two cumulants. These equations involve the density profiles around the TP and 
some tracer-bath cross-correlation functions, whose evolution equations are explicitly 
given in a closed form by resorting to a mean-field-type approximation. In section 
4, we generalize this approximation in order to obtain a closed set of equations for 
the cumulant generating function of the TP position. These equations are valid in 
arbitrary dimension, and they give in principle the whole probability distribution of 
the TP position. As a particular case, we also obtain the equation satisfied by the 
third-cumulant of the distribution. In section 5, we solve the equations obtained under 
the decoupling approximation in the particular case of a one-dimensional lattice, and 
then obtain explicit expressions for the first three cumulants, as well as an implicit 
determination of the cumulant generating function. In section 6, these analytical 
solutions are compared with results from Monte-Carlo numerical simulations that 
exactly sample the master equation. We summarize our results and give an outlook 
in section 7. 
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Figure 1. Model and notations in the two-dimensional case. 



2. Model and master equation 


2.1. Model 


We consider a d-dimensional hypercubic lattice (we denote by a its spacing) in contact 
with a reservoir of particles (see Fig. [^. We adopt a continuous-time description of the 
system. We assume that the particles in the reservoir adsorb onto empty lattice sites 
at a hxed rate f /t*. The particles adsorbed on the lattice desorb back to the reservoir 
with a rate g/r*. The adsorbed particles perform symmetric nearest-neighbor random 
walks, and jump to any of the 2d neighboring sites with a rate l/(2dr*). All the particles 
present on the lattice interact with a hardcore exclusion rule, such that each lattice site 
is occupied by at most one particle. 

We introduce the occupation variable rjr, which takes two values: 1, if the site r 
is occupied by an adsorbed particle, and 0, otherwise. The mean density of the bath 
particles, {pr), is equal to p = //(/ -|- in the long-time limit. However, the number 
of bath particles adsorbed on the lattice is not constant. The case where the number 
of particles on the lattice is conserved can be retrieved by taking the limits / —)■ 0 and 
p —)■ 0 with a hxed value of the density p = //(/ + p) • 

We also introduce a tracer particle (TP). The TP cannot desorb from the lattice 
and it is submitted to an external force, such that it preferentially jumps in the direction 
of the unit vector ei. We suppose that the TP, whose position at time f is a random 
variable denoted by X^, waits a random time distributed accordingly to an exponential 
distribution of mean r, and attempts to jump onto one of the neighboring sites Xt + 

G {e±i, • • ■ , e±d}). For simplicity, we use the notation e_i, = —e^,. The probability 
for the TP to jump in direction z/ is denoted by The equations presented below and 
their solutions are valid for any choice of the jump probabilities. However, it can be 
convenient to assume that the bias originates from an external force F = Fei, so that 
the jump probability in direction z/ writes 


Pu 


e2 


hl3F-e^ 




hl3F-e^ 


(3) 


where (3 = l/(fcB^) is the inverse temperature, and will be taken equal to one. Note 
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that this choice of pv fulhlls the detailed balance condition. After the direction of the 
jnmp has been chosen, the TP hops to the target site if it is vacant; it remains on the 
same site otherwise. 

2.2. Master equation 

We begin by introdncing some anxiliary dehnitions. We dehne rj = {rjr} as the entire 
set of the occnpation variables, which is the conhgnration of the lattice at a given time. 
Next, we dehne the joint probability P{X,ri]t), which is the probability for the TP to 
be at site X at time t with the bath particles in a given conhgnration r]. Next, we dehne 
7]'^’^ as the conhgnration obtained from rj by exchanging the occnpation nnmbers of sites 
r and r + (Kawasaki-type exchange [M] dne to the hop of a bath particle) and as 
the conhgnration obtained from r] by the change pr ^ ^ — Vr (Glanber-type exchange 
[2S] dne to adsorption/desorption events). The time evolntion of the joint probability 
P{X,p;t) is given by the master eqnation : 

d 

2dT*dtP{X, v,t)=Yl E 

fL=l r^X-e^,X 

Odr* 

+ [(1 - Vx) P{X - e^, p; t) - {l - r/x+ej P{X, p; f)] 

+ 2dg Y, [(1 - Vr) P{X, p^-t)- VrP{X, p- t)] 

r^X 

+ 2df Y [VrP{X, 77"; t) - (1 - Pr) PiX, p- t)]. (4) 

r^X 

The hrst term of the right-hand-side of (|^ describes the dihnsion of adsorbed particles, 
the second term corresponds to the dihnsion of the TP, and the third and fonrth terms 
are associated to the desorption and adsorption events of the bath particles. 

If no otherwise specihed, the snm over an index /i rnns over the 2d elements 
{±1, • • ■ , ±d}. In what follows, the brackets (■) denote an average over the TP position 
and bath particles conhgnrations with weight P{X, 77 ; t), and Xt = Xt ■ ei denotes the 
position of the TP along the direction of the external force. 

3. Equations satisfied by the first cumulants 

3.1. Mean position 

The time evolntion of the hrst moment (Xt) of the TP position can be obtained by 
mnltiplying both sides of (|^ by (X ■ Ci) and snmming over all possible conhgnrations 
{X,p). An alternative way to compnte (Xt) is to write that dnring an inhnitesimal 
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time interval At, the TP position Xt evolves according to 


{ Xt + a with probability pi(l - r/Xt+ei) V’ 

Xt-a with probability p_i(l - ?7xt+e_i)v! 

Xt with probability 1 -pi(l - pxt+ef)^ -P-i(l - hXt+e_i)v> 

and take the average of this equation. Both methods result in the following exact 
equation: 

^ ^ {pi [1 - Kr (f)] - p-i [1 - K.t (t)]}, (6) 

where krit) = {rjxt+r) is the probability of having at time t an adsorbed particle at 
position r, dehned in the frame of reference moving with the TP. In other words, kr{t) 
can be thought of as being the density prohle as seen from the moving TP. 

Note that the approximation obtained by replacing all the local average densities 
kr{t) by the global density p in Eq. 0 (that we will refer to as the trivial mean-field 
approximation) yields 

(A',) = -(pi - p_.)(l - p). (7) 

at T 

A more accurate determination of the evolution of the mean position of the TP 
is given by Eq. ([^ and relies on the calculation of the quantities fce±i(^), which are 
the mean density of bath particles at the sites in the vicinity of the TP. This actually 
requires the computation of the density prohle kr{t) for arbitrary r. The evolution 
equations for kr{t) may be obtained by multiplying the master equation (|^ by px+r 
and summing over all the conhgurations of {X,p). We get the following equation : 

2dT*dtkrit) = ^ (V^ - 5r,e^V-ti) kr{t) - 2d{f + g)kr{t) + 2df 

+- 2^Pd ((1 - VXt^ejV^rjXt^r) , ( 8 ) 

u 

where we define the operator acting on any space-dependent function /: 

VJ{r) = f{r + ef,)-f{r). (9) 

Eq. (|^ is not closed with respect to the density prohles kr, but involves the correlation 
functions (?7xt+e^hXt+r)- The evolution equations for these correlation functions can 
be obtained from the master equation (|^, and one can show that they actually involve 
higher-order correlation functions. Consequently, we face the problem of solving an 
inhnite hierarchy of coupled equations for the correlation functions. We then resort to 
an approximation, obtained by writing the occupation variables as pn = (pr) -|- 6pR, 
and by discarding the terms of order (Spr)^. We obtain 

(pXt+rPXt+e^) = {{{PXt+r) + 6pXt+r){{pXt+e^) + 6pXt+e^)) 

^ {PXt+r) (Pxt+e^) 

kr {t) ) 




( 10 ) 

( 11 ) 

( 12 ) 
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which is valid for r ^ 




In the particular case where r = 




Vxt+e^ e {0,1}, we get 

{{VXt+e^Y) = (vxt+ef,) 


recalling that 

( 13 ) 

( 14 ) 


The approximation (12) then relies on the decoupling of the correlation functions 


(rfXt+rPXt+e^)- It can be seen as a mean-field type approximation, and will be gener¬ 
alized later to study the other cumulants of the TP position. It will be referred to as 
the decoupling approximation in what follows. This approximation will be shown to be 
very accurate later on. 


( 16 ) 


Using this approximation in Eq. (|^, we obtain 

2dT*dtkr{t) = Lkr{t) + 2df, 

if r ^ Cj,. For the sites r = with u = {±1, ±2,..., id} we find 

2dT*dtke,{t) = [L + A^{t)]ke,{t) + 2df, 

where L is the operator 

u 

and the coefficients Af^(t) are defined by 


2dr* 

= 1 + - - ie,)*)]- 


( 16 ) 

(17) 

( 18 ) 


The occupation number of the origin is taken equal to zero by convention. Note that 


Eq. (16) represents, from the mathematical point of view, the boundary conditions for 


the general evolution equation (15), imposed on the sites in the immediate vicinity of 


the TP. Eqs. (15) and (16) together with Eq. (18) thus constitute a closed system of 


equations which suffices to compute the density profiles kr{t). These equations were 
first obtained by Benichou et ah, and solved in the case of a one-dimensional lattice |21] 
and of higher-dimensional lattices [271 [26]. 

3.2. Fluctuations of the TP position 

The time evolution of the second moment {Xt'^) is obtained by multiplying the master 
equation by (X ■ Ci)^, and averaging over the TP position and the bath configuration; 
or, alternatively, averaging the balance equation (|^. The details of this calculation are 
given in Appendix A We get 


A(.W)4((X,.ey> 

= y {pi [{Xt) - C/ei (t)] - P-1 [{Xt) - C/e_i (t)] } 

2 

+ — {pi |1 - *:e,(t)l +P-1 [1 - *;«-,(*)]}. 


( 19 ) 
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where grit) = {Xtpxt+r)■ Knowing that 

^(X,)^ = 2 W 




( 20 ) 


and using Eq. (§. we can deduce an expression for the second cumulant of the TP 
position in the first direction: 


^ - {Xtf) = - ^ [PldeAt) - P-lQe.At)] 

2 

+ — {pi [1 - *=•! (<)l + P-1 [1 - *:.-,(*)]} . 


( 21 ) 


where 


grit) = {dXtdpXt+r) : ( 22 ) 

with 5Xt = Xt — {Xt) and 5rjXt+r = Vxt+r — (vxt+r)- It then suffices to determine 
^e±i it) to compute the fluctuations of the TP position. 

In the trivial mean-field approximation, one has kr = p and = 0 for any r. Eq. 
(21) then reduces to 


^ = ^(Pi +P-i)(l -p)- 


r 


(23) 


Note that Eq. (23) corresponds to the standard result for the fluctuations of the position 
of a biased random walker, where the time t is renormalized by the fraction of unoccupied 
sites 1 — p. In the particular case of a one-dimensional lattice where pi +p-i = 1, the 
fluctuations of the TP position do not depend on the bias. This is in contrast with the 
case where the TP has a discrete time evolution [3B]. 

The evolution equations for fjrit) are obtained by multiplying the master equation 
(|^ by dXrfx+r and summing over all the configurations of (X,p). The details of the 
calculation are given in Appendix B We get the following equation 

2dT*dtgrit) = grit) - 2dif + g)grit) 


2dT* 


2dT 


'^p^{5Xtil-rix t+T') 


+ - [Pl ((1 - PXt+e,)gXt +r+ei) — P-1 {i^-Vx t+e_i )VXt +r+e-i )] 


T 

2dT 

T 


-a 


[pi(l - PeJ -P-l(l - Pe_i)] K{t), 


(24) 


We then notice that this evolution equation involves higher-order correlation functions, 
of the form {dXtpxt+rPxt+e^,)■ As previously, their computation would lead to an 
infinite hierarchy of coupled equations. We then write an extension of the decoupling 
approximation (12), obtained by writing pn = -|- dpn and discarding the terms of 

order iSpn)"^. We find 

(dXtPXt+rVXt+e^) ^ (VXt+r) (dXtPXt+e^) + {6XtPXt+r) (VXt+e^) , (25) 
= krit)geAt) + ke At) grit), (26) 
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which is valid for r ^ e^. For r = e^, using the relation (?7xt+e^)^ = Vxt+e,,, we obtain 

= 9.,. (27) 


2dT* 


If r ^ {e±i,..., e±d}, expandind the term of second line of Eq. (24), we obtain 
'^p^{6Xt{l-r]x )^uVx £+r) 


2dT* 


Pm [{^^tVx t+r+e^ ) ~ i^XtPXt+r) — {dXtrjXt+e^VXt+r+e^) + {dXtPXt+ef,VXt+r)] , 

(28) 


and using the definitions of Qr (Eq. (22)) as well as the decoupling approximation (Eq. 

(29) 


(26)), we obtain 
2dT 


T 


Pm - Vx,+e^)Xf,PXt+r) 


2dT* 


T 


[^i‘Sr{t) -9e,(«)V„A^r(«) - ke^{t)V^gdt)] ■ 


( 30 ) 


Then, we gather the first three terms of the r.h.s. of Eq. (24) and obtain 
v9r{i) -‘^d{f + g)gr{t) + {5Xt{l-rix t+e^ )V^r7Xt+r> 


2dT* 


P^9r (P) 


T 

2dT* 

T 


2dT* 


^ n + —— [1 - ke^{t)] ^ Xf^grit) - 2d{f + g)grit) - — '^p^,g^^{t)X^.krit) 


5^PMPe^(t)V/,fc^(£). 


(31) 


The last two terms of Eq. (24) are recast using the decoupling approximation from Eq. 
(12). One obtains 

[Pl ((1 - PXt+ei)PXt+r+ei) “ P-1 ((1 “ PXt+e_i)PX t+r+e-i)] 

2dT* 


2dT* 


-a 


-a 


T 


{pi[l - fcei(£)] -p-l[l - fce_i(£)]} kr{t) 


2dT* 


-a 


[pi(^T’+ei ^ei^r+ei) P—l(^r+e_i ^e_i^r+e_i)] 


2dT* 


-a 


{pi[l - fcei(£)] -P-l[l - fce_i(£)]} krit) 


2dT* 


-a 


{pi[l - fcei(£)]Vlfcr(£) -P-l[l - ke_^{t)]X-ikrif)] . 


(32) 

(33) 

(34) 


Using Eqs. (24), (31) and (34), we finally obtain the following equation for the evolution 
of the correlation functions ^r'- 
~ 2dT* 

2dT*dtgr{t) =Lgr{t) + —^a{pi[l - ke^{t)]Xlkr{t) -p-l[l - ke_^{t)]X-lkr{t)} 

- ^Plk9e^{t)Xf,kr{t), (35) 
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which holds for all r, except for r = {0, e±i,..., e±d}- 

On the other hand, for the special sites r = with v = {±1,..., id}, the term 
in the second line of Eq. (24) can be rewritten as 
2dT* 


T 


- Vx,+e^)^t,VXt+e,) 


2dT* 


T 


2dT* 


T 


Pm (^^4(1 - VXt+e^){VX 
M 

Pm - PXt+eJ(pXt+e,+e^ “ PX^+ej) 


M7^±i/ 


2dT* 

+- Pu - PXi+eJ(pXt+2e, “ PX^+eJ) 

T 

2dT* 

+- P-u {dXt{l - px^+e_j(px^ - PXt+eJ> • 

r 


( 36 ) 

( 37 ) 

( 38 ) 

( 39 ) 

( 40 ) 


The first term can be written in terms of kr and using Eq. (29). The last two terms 
are rewritten using the decoupling approximation (Eq. (26)), the property (pr)^ = Pr, 
and the conventions fco = Po = 0. We obtain 

2dT* 


2dT* 


'^p^ {6Xt{l - PXt+eJV^Px,+e,> 

M 

Pm ['^u9r{t) -g^^{t)Vf^krit) - k^^{t)Vffg^it)] 

2dT* _ _ _ 

i - Pu {92eM - [keAi)92e,it) + ^26. (t)Pe, (t)]} 

T 

2(ir* _ _ _ 

i —^P-y {-9eM + [K.Ak)9eXk) + KXk)9e.Ak)\] 


(41) 


Using Eq. (41), and after straightforward computations, the first three terms of the 
r.h.s. of Eq. (24) yield 


Y - <5e,,e^V_/,) grit) - 2d{f + p)Pe.(t) i Y^l^ (<5Wf(l - PXt+e J V^Px,+e, ) 

M M 

~ 2dT* 2dT* . 

= [L i AAt)]geAk) - Py9eAt)keAt) -X]PMPe,,(t)V/,fce,(t). (42) 


Finally, using again Eq. (34) to rewrite the last two terms of Eq. (24), we obtain the 
equation verihed by ge^{t) for u G {±1, • ■ ■, id}: 

~ 2dr* 

2dT*dtgeAk) =[L i AAt)]g^Ak) + ikeAA -P-i[l - -ikeAk)] 

- PuPeAAkeAA - YPt^SeAA^vkeAt)- ( 43 ) 

r T ^—' 

M 


Eqs. (35) and (43) then form a closed system of equations for the quantities gr{t), 
provided that the quantities ke^{t) are known. Note that Eqs. (35) and (43) are linear 
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in the correlation functions gr{t), which can be written explicitly in terms of the density 
prohles krit). 

The quantities ^e±i (t) can be deduced from these equations, and one can compute 


the evolution of the fluctuations of the TP position using Eq. (21). 


3.3. Stationary values 

We turn to the limit t —>■ oo. We assume that the quantities kr{t) and gr{t) have 
stationary values, so that 

lim dtkr{t) = 0, (44) 

t^OO 

lim dtgr{t) = 0. (45) 


t^OO 


The existence of these stationary values will be shown afterwards. We will use the 
simplified notations: 


kr = lim kr{t), 


t^OO 


gr = lim gr{t), 

t^OO 

= lim A^{t). 

t —^CXD 

We also define the observables: 

V = lim ^ (Xt ), 
t->.oo clt 

1 d 


K= lim ((X^}-(X,f) 
t^oo 2 df ^’ 


( 46 ) 

(47) 

( 48 ) 

( 49 ) 

( 50 ) 


SO that V and K represent respectively the velocity and the diffusion coefficient of the 
TP in the stationary state. Using Eq. the velocity V can be written in terms of the 


functions kr'. 


V = - ke,) - P-1 (l - ke_,)] . 


Similarly, using Eq. (21), the diffusion coefficient can be written 


r = X [pi (1 “ fei) +P-1 (1 - - T (Pi9., -p-i9._,) 


(51) 


(52) 


The stationary values of the density profiles kr (in particular kf,^^ and therefore 


the velocity V) are obtained by solving Eqs. (15) and (16) with dtkrit) = 0. Similarly, 


solving Eqs. (35) and (43) with dtgrit) = 0, one obtains the stationary values of (7e±i 


and the diffusion coefficient K. 

Note that these stationary equations are valid in any dimension, and allow to 
compute the velocity and diffusion coefficient of the TP under the approximations 


presented above (Eqs. (12) and (26)). Their solutions will be presented in the case 


of a one-dimensional system in section 
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4 . 1 . Governing eguations 

In the previous sections, using a decoupling approximation, we were able to determine 
the stationary equations satished by the quantities kr = {rjr) and gr = {dXtdrjXt+r), 
which are involved in the expression of the stationary velocity V (Eq. (51)) and of the 
stationary diffusion coefficient K (Eq. (52)) of the TP. Here, we aim at calculating the 
higher-order cumulants of Xt, defined by 




du’ 


(53) 


u=0 


where the quantity 

T(M;t) = ln(e^“^*> (54) 

is known as the second characteristic function (or cumulant generating function) of 
Using the balance equation ([^, we get the relation 


+ e 


AuXt 


At 


1-(Pl(l - VXt+er) + P-l(l - VXt+e.J 

T 


(55) 


Note that this equation involves two different averages: an average over the direction of 
the step taken by the TP, and an average on the realizations. Eq. (55) leads to 


d _ n (el- - 1) (e»*(l - ^ (e-i- - 1) (e”*)! - px,+=_.)>, 


and, using the definition of T(M;t) from Eq. (54): 
dT 1 d 


dt (e'“^‘) dt 




= ^ - 1 ) 
T ^ 


(einXi) 


+ Pzl (e--- _ 1) 

r ' 


(e“^*r/xt+e i> 


(emXt) 


We define the following correlation functions 


w. 


.(ti;i) = (e'"-*‘px,+r) and wju;t) = Xxf ■ 


( 66 ) 


(67) 


(58) 


(59) 


Finally, we obtain the following evolution equation for the cumulant generating function 
of the TP position: 

^ = n (e- - 1 ) |1 - S.. (m; «)| + (e-1- - 1) [1 - (m; «)] . (60) 
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Assuming that the quantities We^{u]t) and We_i{u;t) reach stationary values when 
t —)■ oo (their existence will be shown a posteriori), the second characteristic function 
has the following asymptotic behavior : 


'l/(M;f) ~ ^{u)t, 

t^OO 


(61) 


with 


$(u) = ^ (e™" - 1) [1 - WeAu)] + ^ - 1) [1 


Pzf 

T 


(62) 


The relation (61) indicates that all the cumulants of Xt are linear in time in the long¬ 
time limit. In particular, this implies that the n-th moment of the rescaled variable 


Zt = {Xt — {Xt))/\J{Xt^') — {Xt)'^ scales as in the long-time limit. All the mo¬ 

ments of Zt of order greater than 2 vanish when f —)■ oo, and Zt is distributed accordingly 
to a Gaussian distribution at large times. 


This calculation allows us to compute the full distribution of Xt in the long-time 
limit. Assuming that the u-dependance of We^^{u) is known, we can derive from the 
previous equations the probability density function (p.d.f.) Pt{x) = Prob[Xt = x] as 
follows. The quantity is defined by 

OO 

^ Pt{x)G^\ (63) 

X= — OQ 


jg the Fourier transform of the p.d.f. Pt{x), which can be obtained by the 
inverse Fourier transform: 


Pt{x) 



du 




(64) 


and, using the expression of the long-time limit of \h(f) (Eqs. (61) and (62)), 


Pt{x) ~ [ ^exp | —(e‘“‘^ - 1)[1 - Wei(M)] + ^^(e - 1)[1 - We_i(M)] - mxl. 

ZTT I T T J 

(65) 

Consequently, it suffices to determine the u-dependence of uje±^(u;t) to obtain the 
p.d.f. of the TP position. In what follows, we establish the evolution equations for 
the quantities Wr{u;t) starting again from the master equation (|^. 


4-2. Evolution equations of the quantities Wr{u;t) 


The evolution equation of the correlation functions Wr{u]t) (defined by Eq. (59)) 
can be obtained by multiplying both sides of the master equation (|^ by the quantity 
and averaging with respect to the bath configuration r] and the TP position 
Xt- Extending the method used to derive the evolution equations of the correlation 
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functions gr{t) starting from the master equation Q, it is found that Wr{u] t) obeys the 
following exact equation (see Appendix B): 

2dT*dtWr{u] t) = Wr{u; t) - 2d{f + g)wr{u] t) + 2df (e‘“^*) 


2dT* 


r 


2dT* 


T 


- Vx,+e,)Vf,r]x,+r) 

Pe (e'“'" - l) (e™^*(l - gxt+e,)gxt+v+^,) 


( 66 ) 


e=±l 


We then make the following decoupling hypothesis: 

((e“- - 'te.+,<?x,+e„> =: ((e“* - (e'”*» (i)x,+e,> 

+ ((e‘“* - (e»*)) (Vx,+r ). 

which is valid for r ^ e^. This is equivalent to 

(D^^^Pxt+rVxt+e^) = Wr{u;t)K^{t) + kr{t)w^^{u; t) - (e™^‘) K{t)K^{t). 

For r = e^, one gets 

(e'”*(r;x,+e„)^> = (e'“*r;x,+e„) 

= «>.,(«;()■ 


( 67 ) 


( 68 ) 


( 69 ) 

( 70 ) 


This decoupling approximation is an extension of the approximations (12) and (26): it 
is obtained by writing the occupation variables as pu = {pr) + Spu, and discarding the 
terms of order {6pit)^. Note that expanding Eq. (68) at order 0 and 1 in u, we retrieve 
the decoupling approximations made for the correlation functions (pxt+rPXt+e^) (Eq- 
(12)) and {dXtpxt+rVxt+e^) (Eq. (26)). Using this approximation in Eq. (66), and 
following a procedure similar to the one used to derive Eqs. (35) and (43) from Eq. (24) 
we obtain for r ^ e^, 


2dT*dtWr{u]t) = LWr{u]t) + 2df (e'“^‘) + Yp^^ [^e^{t) (e‘“^‘) - We^{u]t)\ V^kr{t) 

H- Y {'^r+e, [u] t) [1 - ke, (t)] + /Cr+e, {t) ) k^, {t) - [u] f)] } . 

^ t=±l 

(71) 

For r = Ey, the evolution equation becomes : 

2dT*dtWeAu]t) = [L + A„{t)]we,{u]t) + 2df (e™^*) 

OrJr* _ 2(7t* 

+- YPf^ -We^{u-,t)] Vf,ke,{t) + - p^ [KAt) - We,(MU)] KAt) 


2dT* 


Y (e™"'" - l) {We.+e,(M;t)[l - k^AA] + ke.+eAA k^AA “ ^eA^A)] } 


e=±l 


( 72 ) 
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These equations are conveniently written in terms of the variable Wr{u]t) by noticing 
that : 


dwr 


1 dWr{u]t) Wr{u]t) d 


dt (Q^uXt ^ Qi 

1 dwr- _ / 

Wr{U]t) — 


fg\uXt) Qf 


(^QiuXt ^ Qf 

1 dWr 

{e^uXt ^ 


dt 


- w^(u;t) - l) [1 - (e - 1) [1 - Wei(u;t)] 


(73) 


Finally, we divide the evolution equations (71) and (72) by and obtain the 

evolution equations for Wr(u;t), for r ^ : 

2dT*dtWr(u;t) = LWr(u;t) + 2df + ~ ^e^(u;t)] Vf,kr(t) 


2dT* 


Pe - l) {VeWr{u]t) - fee, (t) [Wr+e, (m; “ 4+e,(t)] “ [/^r+e, (^) “ 


e=±l 


and for r = : 


2dT*dtWf,^{u] t) = [L + A^{t)]we,{u] t) + 2df 
2dT* . 2dT* 

+ — -We^{u]t)\ V^ke,{t) + —^Pu [ke,{t) -We,{u;t)]ke^{t) 

2dT* x 

+- Y1 “ A;e.+e,(t)] - Wg, (w; f) [fce,+e 

^ e=±l 


M 


These equations can in principle be solved in the stationary limit t —>■ cx), by setting 


dtWr{u]t) = 0 in Eqs. (74) and (75) and obtaining the values of Wr{u) satisfying these 
stationary equations. In particular, this allows us to obtain the u-dependence of the 
functions We±i{u) and to deduce Pt{x) from Eq. (65). These equations are valid in any 


dimension, and their solution gives the cumulant generating function and the probability 
distribution function of the TP position. 

This resolution will be made explicit in the case of a one-dimensional lattice in 


section 6.5 We also notice that the functions Wr{u]t) can be expanded in powers of 
u to compute higher-order cumulants. In particular, we give in the next section the 
evolution equation satisfied by the third cumulant of the position of the TP. 


f.3. Third-order cumulant 

In this section, we study the third cumulant of the distribution of Xt, which characterizes 


its skewness. We use Eqs. (74) and (75), describing the evolution of Wr {u',t), to calculate 


(«;*)]}. 

( 74 ) 


( 75 ) 
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the third order cumulant. We dehne the coefficient 7 by the relation : 

t^oo 6 dt' ' 


( 76 ) 


Recalling the dehnition of from Eq. (54), we get the following expansion in powers 
of m: ^ ^ 

*(«) in {X,) + {{X, - {X,)f) + {{X, - {X,)f) + ... (77) 


6 


In the long-time limit, and using the dehnitions of V (Eq. (49)), K (Eq. (52)) and 7 

( 78 ) 


(Eq. (76)), one gets 


dd' 


lim ^ = iuV + (m)^R' -|- (m)^7 -|- . 

t—>oo q/) 


We dehne the correlation function mr{t) by the relation 

fhrit) ^ ((W - {Xt)Yvx,+r) - K{t) - {X,f] , 

SO that the expansion of Wr(t) in powers of u writes 


( lU ) 

Wr{t) = kr{t) + iugr{t) + X- — fhr{t) + O [u^) , 


(79) 


(80) 


where we used the dehnitions of krit) = {rjxt+r) and gr{t) = {5Xt5rjXt+r)■ Expanding 
both sides of Eq. (60) up to order 3 in m in the limit t ^ 00 and using Eq. (78), one 
gets the following expression for 7 : 

7 = 


PlCr 

^a^{l - fcei) - 


1 

to 1 

P-iO- 

[1 9. 

^e_i) T 2^9e-i 

1^ 

-mp , 

2 

r 

T 








(8 


According to Eq. (80), the general evolution equations for Wr{u]t) (Eqs. (p^ and 


(75)) expanded at order 2 in m then give the evolution equations of mrit). We get: 
• for r ^ Gi, ■. 


2dT*dtmr{t) = Lmr{t) 


2dT* 


T 


^kr{t) 


H-p,(T{2e[(l - keXt))X elfrit) - ekr{t)] + Ct( 1 - fee, (t) ) (t) } • 

( 82 ) 


e=±l 


for r = Gy ■. 

^ ~ ^ 2dT* , _ 2dT* ^ 

2dT*dtm^^{t) = [L + A^(f)]me,(t)-^ f^k^^jt) - —p^m^^{t)ke,{f) 


T 


r 


2dT 


h- 5^ p,cr{2e[(l - keM)Xe9eAt) -PeM^ekeAt)] + ^^(1 “ (t)) V,/Ce. (t) } • 


e=±l 


(83) 


In the stationary limit, one computes from Eqs. (82) and (83) the quantities m 


e±i 


and then the coefficient 7 from Eq. (81). This will be made explicit in the case of a 
one-dimensional lattice in section 15.1 
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5. First cumulants and distribution of the TP position in one dimension 


In this section, we focus on the one-dimensional version of the general model presented 
in section |2.1[ This situation is related to a number of lattice models of interacting 


particles which have been widely studied in the mathematical and physical literature. In 
particular, in the situation where the number of particles on the lattice is conserved and 
where the TP is not biased, the model corresponds to the well-known single file problem, 
for which several results have been derived exactly [STiiaHiEHiiioiiiiiiiaiiaiiii]. 

In what follows, we consider the equations derived in the previous section in the 
particular case of a one-dimensional lattice. We first recall the solutions of the equations 
satisfied by the density profiles kr that were obtained in previous studies, and give a 
detailed resolution of the equations verified by the correlation functions Qr and Wr- 


5.1. Solution of the equation on kr in one dimension 

The solutions of the equations verified by kr have already been presented [21], and we 
recall them here for completeness. In one dimension, we adopt the simplified notation 
kr = knei = kn- The stationary limit of the general equation verified by the density 
profiles (Eq. (HI) is then a second order recurrence relation on the quantities kn- Its 
solution has the following form: 


kn 


p + K+ri^ for n > 0, 
p + K^rf^ for n < 0, 


(84) 


where 


n = 
2 


Ai + A-i + 2 (/ + g) T \j -|- A-i + 2 {f + g)^ — AAiA^i 


2 A. 


while the amplitudes K± are given respectively by 

Ai — A i 


and 


K^ = p 


K_ = p 


dl-i-Ain’ 
Ai - kl i 


(85) 


( 86 ) 


(87) 


A-i/r2 — Ai 

We notice that r2 > ri, which indicates that the density profile behind the TP decreases 
slower than in front of it. We also note that iP+ > 0 and < 0 which indicates 
that there is a jammed region ahead of the TP, and, on the contrary, a depleted region 
behind it. 

One obtains a closed sed of two non-linear equations determining implicitly the 
parameters Ai and A_i, from which one can compute the TP stationary velocity, related 
to A±i through 

( 88 ) 


a 
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Substituting Eq. (84) into the definition of (Eq. (18)), we find 


Ai — 1 + 


A_i — 1 + 


2piT* 

T 

2p-ir 

T 


P- P 


p- p- 


Ai-A_i 
A_i/ri - AiJ ’ 
Ai - 2l i 


(89) 

(90) 


A_i - Air2_ 

For a given set of parameters (/, g, a, r, r* and pi), the numerical resolution of this 


system leads to the values of Ai and A_i and then, using Eq. (88), to the value of the 


stationary velocity of the TP. This approximated value of the stationary velocity will 


be compared to numerical simulations in section 6.2 


5.2. Solution of the eguation on pr in one dimension 

We now go one step further and determine the diffusion coefficient K. This in turn 
requires the knowledge of the functions fjr- For simplicity, we adopt the notations 


Qr = 5'nei = On- Usiug the expressions of kn from Eq. (84), the general equations 


satisfied by Pn (Eqs. (35) and (43)) become: 
• for n > 1: 


for n < —1: 


■^i{.9n+i 9n} T A_i(^gn_i Qyfj 2(y T 9')9n 

+ ^cr l^piK+rf - p - K+n - (ri - 1) 

- P-iK+rf - p - K+ri + (rf ^ -1)1=0, 


Al{9n+1 9n) T A—ii^Pn—l 9n) 2(/ + p)pn 

+ |piA’_rJ (^1 - p - A’+ri - (r2 - 1) 

- p-i/f-r; (l - p - K+ri + ^) (rp - 1)1 = 0, 


(91) 


(92) 


Pi and p_i may be computed using the boundary conditions given by Eq. (43): 


^92 - 9i ( ^-1 + 2(/ + 9 ) + ‘^PiiP + K+rl)j + ^p_i(p + K+ri)g_i 


2 D 


api{l - p - K+ri){p + K+rl) 


2t* 


+-c^(Pi(l-p-^+n) -p-i{l- p- K_r^ ))(p + iP+ri), 

T 

p4i + 2(/ + p) + ^p_i(p + K_rf^)j + ^Pi(p + K_rf^)gi 
= ^ap_i(l - p - K_rf^){p + K_rf^) 

+—a(pi(l - p - iP+ri) - p_i(l - p - K_rf^)){p + K_rf^). 

T 


(93) 


( 94 ) 
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The general solution of Eqs. (91) and (92) can be written: 


and 


g-n = ar^ - 


gn = /3r^ - 


W 


Airi - A_ir^ 

W 


^nr"^ for n > 0, 


zi^nrlf for n < 0, 


Air2 - v4_ir2 

where a and ft are constants to be determined, and where 


W = K+ - a <( Pi 


2t* 

W = K_ -a 

r 


Pi 


1 ir 3i 

I- p- K+ri - 

a 


1 ir 

I- p- K+ri - 

a 


(n - 1) -p_i 


1 ^ 3-1 

1 -p -+- 

ri a 


(r/ 


(r2 - 1) - p-i ( 1 - p - — + — 

ri a 


(95) 


(96) 


- 1 ) 

(97) 


(’’2 - 1) 

(98) 


Substituting Eq. (95) into Eq. (93), Eq. (96) into Eq. (94), and writing Eq. (95) for 


n = 1 and Eq. (96) for n = —1, we obtain a linear system of four equations satisfied by 


a, ft, gi and g-i, which is straightforward to solve. The explicit expressions of gi and 
p_i are given in Appendix C Note that they rely on the determination of the quantities 


K± and ri,r2, which are determined numerically for a given set of parameters with the 


method detailed in section 5.1 Finally, for a given set of parameters, one can deduce 


the values of g±i and the value of the diffusion coefficient using Eq. (52). Note that this 


calculation also gives access to the spatial dependence of the cross-correlations functions 


gn through Eqs. (95) and (96) 


In section 6^, we investigate the dependence of K on the different parameters of the 
problem, and we confront the analytical prediction from the decoupling approximation 
to results from numerical simulations. We give a first insight into the understanding 
of the counter-intuitive non-monotonic dependence of K over the density p that was 
described in [31]. 


5.3. Solution of the eguation on rur in one dimension 

We finally solve the equations satisfied by the correlation functions me±u from which 
we will compute the coefficient 7, related to the third cumulant of the distribution, and 


defined by Eq. (79). Starting from the general equations verihed by and valid in any 


dimension (Eqs. (82) and (83)), we study the one-dimensional case. For simplicity, we 


write ror = ronei — ^n- The quantities are the solutions of the equations presented 


below: 


for r ^ e^, using Eq. (82), one gets : 


- nrtn) + A_i{mn-i - rUn) - 2(/ -F g)mn = S{n) 


(99) 
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where S{n) can be expressed explicitly in terms of the fnnctions kn and gn, 
determined respectively in sections 5.1 and 5.2, For n > 0, we write S{n) nnder 
the following form : 

^(n) = (C+ + nCf)r^. (100) 


solntion of Eq. (99) reads: 


The explicit expressions of Cf and Cf are given in Appendix D Conseqnently, the 

fhn = Trf + (a+n^ + b+n)rf, (101) 


with 


2 AiVi — A_ir^ ^ ’ 

= -- — - - a+(Airi + A_irf^]. 


Airi - A_ir{ 

For n < 0, a similar resolntion leads to: 


with 


1 

a_ = - 


uin = Srlf + (a_n^ + 

Cf 


b_ = 


2Air2 - A_ir2 

- a_(Air2 + A_irf^], 


( 102 ) 

(103) 

(104) 

(105) 

(106) 


and where the explicit expressions of and C 2 are given in Appendix D 


for r = Cl, we obtain from Eq. (83): 
0 = Aim2 — fhi 

with 


2t* 

T 


2 r 

T 


(fl = 


2 dD 


T 


2 dD 


T 


-P- 


for r = e_i we obtain from Eq. (83): 


2t\ 

T 


+ cpi, 

(107) 

p2 - ki){l - ki)} 


a^ki{l - k-i)} . 

(108) 

2 t* 


- Pimik_i + 

(109) 


with 


P-i = 
2 dT 


2 dD 


-P 


'_i {-2cr [-fi'-i(l - ki) + gik-i] - a'^k-i{l - ki)} 


ZiCLT 

^- P-l {2cr [(5-2 - 5-l)(l - k-i) -g-i{k^2 - fc-l)] + - ^-1)} - 


( 110 ) 
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for 77, = —1, and considering the bonndary 
one obtains a linear system of fonr eqnations 

with unknowns mi, m_i, F and 5 : 


Writing Eq. (101) for n = 1, Eq. (104 
conditions given by Eqs. (107) and (109), 


^ T\ 0 ^^13 


/ r \ 


(yA 

0 r2 ^ -^23 M2A 


5 


Y2 

Air\ 0 M33 M34 


rhi 


Ys 

Y 0 A—\t 2 ^ -/W43 A/44y 


\fh-i) 


vJ 


( 111 ) 


where the expressions of the quantities Mij and Yj are given in Appendix D Finally: 


(4ei M ^i2M24){MiiY2, + M^iYi) 


-|-(M34Mii — M^iMi4)[M42Y2 + M22F4)] , 

(112) 

tri-i Y[42M2‘i){MiiY2, + M^iYi) 


— Mi^M^i){M42Y2 + A722F4)] • 

(113) 


The procedure to compute the coefficient 7 for a given set of parameters is the 
following. With the method presented in section 5.1, one can compute numerically the 
quantities iF±, ri, r2 and k±i for a given set of parameters. Using the analytical expres¬ 
sions of ^±1, a and (3 in terms of these quantities given in Appendix C, one computes 
fh±i with Eqs. (112) and (113). The coefficient 7 is deduced from its dehnition (81). 


5.4- Solution of the equations on Wr in one dimension 


We now turn to the resolution of the equations satished by the correlation functions Wr 
in the specihc case of a ID lattice and in the stationary limit. 


Starting from Eqs. (74) and (75) and assuming that there exists non-trivial 


stationary solutions, one gets the following equations satished by Wn = Wnei ■ 

BiWn+i - B 2 Wn + B^Wn-i = -{B^K+r"^ + B^) for n > 1 (114) 

BiWn+i - B 2 Wn + B^Wn-i = -{CiK_rl4 + B 5 ) for n < -1 (115) 

D 1 W 2 — D 2 W 1 + D 3 W -1 + D 4 + D^wl -f DqWiW-i =0 for n = 1 (US) 

EiWi — E 2 W -1 + E 3 W -2 + E 4 + D^wfi + D^wiW-i =0 for n = —1 (117) 

where the expressions of the different coefficients Bi, Ci, Di, Ei are given in Appendix 
1^ Eqs. (114) and (115) are associated to the characteristic equation 


BiX^ - B2X + 53 = 0 , 


(118) 


<?1,2 = 


B2 ± ^Bl - AB^B^ 

2 Bi 


which has the solutions 


(119) 
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qi = ri + 0{u) and q 2 = r 2 + 0{u). 


( 120 ) 


At order zero in u, Eqs. (114) and (115) are equivalent to the stationary limit of Eqs. 
(15) and (16), so that their solutions must coincide at this order. Consequently, the 

( 121 ) 


general solution of Eq. (114) is of the form 


Wn = a+gi 


The particular solution of Eq. (114) is easily calculated, and for n > 0, we find 

B,K+ 


Wn = a+gi - 


-1 ' 1 


r?- 


BiVi — B 2 + B^r^ 

With similar arguments, we find for n < 0 

C 4 K. 

Brr2-B2 + B^rf^ ^ 


B, 

Bi — B 2 + B 3 


( 122 ) 


B, 

Bi — B 2 + S 3 


(123) 


Finally, writing Eq. (122) (resp. Eq. (123)) for n = 1 (resp. n = —1) and using 
Eqs. (116) and (117), we find the following nonlinear system of four equations whose 
unknowns are Wi, W-i, «+ and «_ : 

B,K+ 


S, 


wi = a+qi - 


tc_i = Q!_g 2 
Si 


-1 


BiVi - B2 + B^rl 
CaK. 


-1 


n - 


Bi — B 2 + S 3 

S 5 


a+qf 


Sir2 - S 2 + B^r^ 
BaK+ 2 

—n 


- 1'2 


Vo ^ - 


Si- 

S 5 


Bo 


S 3 


Siri - S2 + B^rf^ " Si - S2 + S3 
S2u;i + D3W-1 + S4 + D^wf + DqWiW-i = 0 
-2 CaK_ _2 

CK-On-l-rn 

Sir2 — S 2 + B^r2 


(124) 


EiWi — E2W-1 + S3 


S. 


Si — S 2 + S 3 


+ S 4 + Dewfi + D 5 W 1 W -1 = 0. 


The numerical resolution of this system of equations for specific values of u allows us to 
calculate wi and W-i as functions of m, and to deduce the stationary cumulant gener¬ 
ating function T as a function of u. Using Eq. (65), one can calculate the probability 
distribution Pt{x), valid in the asymptotic regime t — )■ 00 . 


In the next section, we analyze the solutions obtained for the velocity V, the 
diffusion coefficient K and the coefficient 7 to study their dependence on the different 
parameters of the problem. These results, which were obtained using our decoupling 
approximation, are confronted to numerical simulations which exactly sample the master 
equation (|^. 
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6.1. Algorithm and numerical methods 


In order to verify the accuracy of the approximation involved in the computation of 
the cumulants of the TP position, we perform numerical simulations. We use a kinetic 
Monte-Carlo (or Gillespie) algorithm [l5l H6] in order to get an exact sampling of the 
master equation (|^ describing the dynamics of the system. The details of the numerical 
methods are given in Appendix F 


6.2. Velocity 

For completeness, we present results for the velocity of the TP, which had already been 
presented in [23] . We study here the terminal velocity reached by the TP as a function 
of the density p, for different values of the bias. As p is in fact hxed by the values of / 
and p, we decide to vary / for different values of g in order to explore the whole range 
of parameters. Results are presented in Fig. As expected, the velocity of the TP is a 
decreasing function of the bath density. We also confront the result from our decoupling 
approximation to the trivial mean-field solution (dashed line). 

The discrepancy between the results from numerical simulations, which correspond 
to an exact sampling of the master equation (|^, and the solution obtained using the 


decoupling approximation (12) is very small, and the agreement is particularly good 


close to p = 0 and p = 1. The decoupling approximation is then very accurate for the 
estimation of the velocity of the TP. 


6.3. Diffusion coefficient 

6.3.1. Results In a similar way, we study the diffusion coefficient as a function of the 
density p, for different values of pi and g. We confront the analytical predictions to 
results from numerical simulations in Fig. The dependence of the diffusion coefficient 
on the density of particles was hrst investigated in |3T|. In the simulations results 
presented in that publication, the numerical errors were underestimated, and we present 
here refined results. 

The analytical predictions as well as the numerical simulations reveal the existence 
of a striking effect: the diffusion coefficient can be a non-monotonic function of 
the density of bath particles. Conterintuitively, increasing the density of particles 
surrounding the TP may increase the TP diffusion coefficient in some range of 
parameters. This result is surprising, since one naturally expects that the diffusion 
coefficient will be maximal when there are no hardcore bath particles (i.e., when p = 0). 
Consequently, this means that the diffusion of the TP may be enhanced by the presence 
of bath particles on the lattice. Moreover, we observe that for given values of pi and 
g, this extremum of the function K{p) may appear if the bias pi — p_i is large enough. 
This effect could be investigated in experimental situations (e.g. in microrheology) and 
could have interesting applications. 
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Figure 2. Stationary velocity V of the TP as a function of the density for different 
values of the desorption rate g obtained from numerical simulations (symbols) and 
from the decoupling approximation (lines). The bias is pi — p-i = 0.96 (top) and 
Pi — p_i = 0.6 (bottom), the waiting times are r = t* = 1. The dashed line is the 
trivial mean-field solution V = ^(pi — p_i)(l — p). 


The discrepancy between the resnlts from nnmerical simnlations, which correspond 
to an exact sampling of the master eqnation (|^, and the solntion obtained nsing 
the deconpling approximations (Eqs. (12) and (26)) is small, except in the domain 
0-05 ^ p ^ 0.15. However, the approximate solntion still gives a good qnalitative de¬ 
scription of the evolntion of K. The agreement is particnlarly good close to p = 0 and 

p= 1. 
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Figure 3. Stationary diffusion coefficient K of the TP as a function of the density for 
different values of the desorption rate g obtained from numerical simulations (symbols) 
and from the decoupling approximation (lines). The bias is pi — P-i = 0.96 (top) and 
Pi — p_i = 0.6 (bottom), the waiting times are r = t* = 1. The dashed line is the 
trivial mean-field solution K = |^(1 — p). 


In what follows, we find the criterion on the parameters g and pi allowing the 
emergence of a maximum for K{p). In addition, we show that the non-monotonicity 
of the diffusion coefficient is actually correlated to the non-monotonicity of the cross¬ 
correlation functions in the domain n < 0 (i.e. behind the TP). 

6.3.2. Criterion for the existence of a maximum value of K In this section, we 
determine an explicit criterion for the existence of this maximum. This is equivalent to 
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determine the parameters for which K{p) has a positive derivative at the origin. We 
will then solve the equation: 

^ =0. (125) 

dp p=0 

If g is fixed, as the density p is related to / and g through p = f/{f + g), this is 
equivalent to considering K as a function of / and solving 


dK 


= 0 . 


(126) 


1=0 


In what follows, we obtain the leading order term of K in an expansion in powers of 
/, the other parameters being constant. For simplicity, we introduce the quantities 
t' = T*fr and 6 = pi — p_i. Assuming that the quantities A±i have the following 
expansions 

2l±i = 21™ + A'Pj + 0{f), (127) 


and using Eqs. (89) and (90), we obtain 


A±i = l+(l±(5)r' 

f^O 


1±(5 


+ g(l + t) ± + g^ + 2g(l + t')I f+0(f) 

(128) 


We can deduce from the expansions of A±i the expansions of K±, ri and r 2 using Eqs 


, (86) and rt87|. We then obtain the expansions of k±i and g±i in powers of / using 


the results from sections 5.1 and 5.2, and finally an expansion of K with Eq. (52). 


These general expressions are too lengthy to be reproduced here, but we give them in 
the case r' = 1, which is the case we considered in our simulations: 


K = — 
f^o 2r 




(129) 


where 


N{g, 6) = 2 [{6^ - 4 )/ + {36^ - 16 )/ + 6^{36^ -56^ + 4)g + 6^6^ - 36^ + 4 )] 

-V^^ + 9^ + ^9 [ 2 ( 5 " - 4)(52 - 3 )/ + 52(^4 _ ^2 ^ ^ 4^4(^2 _ 2 )] ^ ( 130 ) 

Dig, 6) = 2(5 - 2 ) 2(5 + 2 )^( 5 ^ + / + 4 //. ( 131 ) 

For any value of 5 and g, Dig, 5) > 0. For a given value of the bias 5, the critical value 
of g (denoted by g^f) allowing K to reach a maximal value is then the solution of the 
equation: 

7V(^c,5) = 0, (132) 


which can be determined numerically. We present the numerical solutions of Eq. (132) 
in Fig. 1^ (curves in blue). On this figure, we also give the solutions obtained for other 
values of r'. We conclude that, for a fixed value of the bias, K is nonmonotonic if the 
desorption rate g is large enough (or, for a fixed value of the desorption rate, if the bias 
is large enough). 
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0 0.2 0.4 0.6 0.8 1 

S 

Figure 4. Critical value of the desorption rate gc as a function of <5, for different values 
of the ratio t' = T*/r, with cr = 1 and r = 1 obtained by the study of the behavior of 
K at p ^ 0 (curves in blue). The region below the curves correspond to the range of 
parameters where K{p) is nonmonotonic. The critical value of the desorption rate 
for which the cross-correlations 5 have the property g-i > g -2 is represented in green. 
The region below the curves correspond to the range of parameters where g_i > 5 - 2 ■ 


6.3.3. Influence of the cross-correlations functions on the non-monotonicity We now 
aim to give a hrst insight into a better understanding of the non-monotonicity of the 
diffusion coefficient. We first notice that the diffusion coefficient K (Eq. (52)) may be 
separated into three contributions: 


K = Kmy + K^ + K 2 , 

(133) 

with 




0 ^ 

-^MF = -(1 ~ P)) 

T 

(134) 


Ki = - [pflki-p)+p-i{k-i-p)f 

T 

(135) 


K 2 = - —{pigi-p-ig-i). 

T 

(136) 


The first term is the trivial mean-held approximation of the problem, obtained by 
taking the average local densities kr equal to p and the cross-correlation functions = 
{SXtdrjXt+r) equal to zero. Ki may be seen as a contribution from the inhomogeneous 
density prohles, and K 2 a contribution from the cross-correlations 'g±i. For a given set of 
parameters which gives rise to a non-monotonic behavior (pi = 0.98 and g = 0.15), we 
plot in Fig. as well as the three contributions. The origin of the non-monotonicity 
of K with respect to the density p can then be attributed to K 2 - 

The cross-correlation functions Ijn have nontrivial behaviors with respect to the 
distance n to the TP. In particular, they appear to have a non-monotonous behavior in 
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Figure 5. Contributions in the expressions of if as a function of the density, for the 
parameters pi = 0.98 and g — 0.15. 




Figure 6. Cross-correlations functions g as a function of the distance to the tracer n 
for g = 0.2 (left) and g = 0.6 (right), and for different values of the density p. In both 
cases. Pi = 0.98. 


the domain n < 0, i.e. behind the TP, for some values of the parameters. On Fig. 
we plot the functions for different values of the density p, for two sets of parameters: 
one for which the diffusion K is known to be a non-monotonic function of the density 
(pi = 0.98 and g = 0.2), and one for which it is monotonic (pi = 0.98 and g = 0.6). On 
Fig. 0 we plot the functions gn for different values of the desorption parameter p, for 
p = 0.01, for Pi = 0.98. The non-monotonicity of with the distance to the TP then 
seems to be correlated with that of K with the density, as it occurs for g small enough. 

In order to get a more quantitative comparison of these two effects (the non¬ 
monotonicity of K with respect to the density p and that of with respect to the 
distance to the TP n in the domain u < 0), we hrst determine, for a hxed value of the 
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Figure 7. Cross-correlations functions g as a function of the distance to the tracer n 
for p = 0.01, Pi = 0.98, and for different values of the desorption parameter g. 


bias Pi — P-i, and at leading order in /, the critical valne of giving a non-monotonic 
behavior for Qn, i.e. for which ^_i > ^_ 2 . 

Using the small / expansion of A±i (Eq. (128)), we easily dednce from the 
dehnitions of Qn in the domain n < 0 (Eq. (96)) an expression for g^i — ^_2 at leading 
order in /. The general expression is too lengthy to be given here. We present it in the 
particular case where r' = 1: 




(137) 


with 


Nfg, 5) ={5 - g) [(1 + 5)g^ - 2{5^ -25- 2)g + 5^5- 1)] 

- ^5^ + g^ + Ag [(1 + 5)g^ - (35 + 2)(5 - l)g + 5^(5 - 1)] , (138) 

D\g, 5) =2g^5^ + g^ + Ag[ [2g^ + (25 + A)g + <5(5 + 2)] ^5^ + g^ + Ag 

+ 2g^ + 2(5 + 4)/ + 2(5 + 2){5 + l)g + 5^{5 + 2)}. (139) 

Then, for a hxed value of the bias 5, the critical value of g canceling fj-i — fj -2 is the 
solution of the equation 

7V'(^',5) = 0. (140) 

The numerical solutions for g'^ as a function of <5, for different values of the ratio r' are 
represented in Fig. and confronted to the values gc obtained by the criterion on K. 
The two functions gc and g'^ are comparable as long as r' is not too large. We also show 
that their expansions for 5 —)■ 0 are identical at leading order, and that they go to zero 
as (vTf -|- l)5^/8. 

This study indicates that the emergence of a maximal value for K{p) for p > 0 is 
correlated to the observation of a minimal value in the cross-correlation functions in 
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In [Appendix G[ for given values of the jump probability pi and of the desorption 
parameter g, we study the range of density [ 0 , p\ for which the diffusion coefficient is 
greater than 1/2, which is its value when there is no bath particle. We also deter¬ 
mine the range of density [Ojp'] for which — ^_2 < 0. The range of parameters for 
which K > 1/2 is then shown to be correlated to the range of parameters for which 
g-i - 9-2 < 0. 


These two studies (from sections 6.3.3 and Appendix G) show that the non¬ 


monotonicity of the diffusion coefficient with respect to the density (emergence of a 
maximum value of K for p > 0) and the non-monotonicity of the cross-correlation 
function behind the TP are correlated. A more detailed study of the cross-correlation 
functions whose behavior affects the fluctuations of the TP position, could allow us 
to have a more physical understanding of the phenomenon highlighted in this section. 


6.4- Third cumulant 

The third cumulant of the distribution of Xt gives information about its asymmetry. 
We introduced earlier the coefficient 7 , dehned by 


7 = lim {{Xt 

t^oo 6 dt 


{x.Yf) ■ 


(141) 


With this dehnition, if 7 > 0 (resp. 7 < 0), the distribution of Xt is expected to be 
skewed to the right (resp. to the left). Let us recall the situation of a biased random 
walker on a lattice in the absence of exclusion interactions. If the particle is more 
likely to jump to the right {pi > p-i), its third cumulant will be positive, indicating a 
distribution skewed to the right. Here, we study the influence of the presence of bath 
particles and of the different parameters of our model on the sign of the third cumulant 
ofW. 


We use the solutions for m±i (Eqs. (112) and (113)) obtained from the decoupling 


approximation in order to compute the coefficient 7 from its dehnition (Eq. (81)). For 
two values of the bias {pi = 0.98 and 0.8), we study the coefficient 7 as a function of the 
density p, for different values of the desorption rate g. The solutions of the equations 
obtained from the decoupling approximation are compared to results from numerical 
simulations. The curves are presented in Fig. For high values of the desorption rate, 
7 (p) is monotonic and decreases when the density of bath particles increases. However, 
for small values of the desorption rate, the function becomes nonmonotonic, and one 
observes the emergence of a minimum and a maximum value, different from the trivial 
extrema at p = 0 and p = 1. Finally, if g is small enough, there exists an interval of 
density for which 7 becomes negative, which means that the distribution of Xt may 
actually be negatively skewed, in opposition with the situation where there is no bath 
particles on the lattice. 
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Figure 8. Coefficient 7 defined in Eq. ( 141 ) as a function of the density for different 
values of the desorption rate g obtained from the decoupling approximation. The bias 
ispi—p-i = 0.96 (top) andpi—p_i = 0.6 (bottom). The waiting times are r = r* = 1. 
The dashed lines are the trivial mean-field solutions: fy(_Pi — p_i)(l — p). 


Confronting these predictions with results from numerical simulations, we show 
that the decoupling approximation offers a good prediction of the behavior of the third 
cumulant of the TP position in a wide range of parameters. 


6.5. Cumulant generating function and propagator 

For a given set of parameters, we solve the system of equations which determines w±i 
implicitly (Eq. ( |124 )), and we obtain the numerical values of w±i{u) for u varying in 
[—TT, vr]. We represent the real and imaginary parts of w±i{u) as functions of u in Fig. 
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Figure 9. Real and imaginary parts of the generalized correlation functions w±i{u), 


obtained from the resolution of the system (124), for the parameters p = 0.9, g = 0.15, 
Pi = 0.98, T = T* = cr = 1. 


1^ The p.d.f. Pt{x) in the long-time limit is obtained from Eq. (65). 

We compare the results from this calculation with data obtained from Monte Carlo 
simulations for a given set of parameters (see Fig. 10). We observe a good agreement 


between the analytical prediction obtained from the decoupling approximation and the 
results from numerical simulations. We see that the prediction from the decoupling 
approximation tends to be shifted to the right for large times: this is expected from 
the analysis of the velocity of the TP, which was shown to be overestimated by the 
approximation (Fig. [^. 


As emphasized in section 4.1, the rescaled variable Zt = (Xj — (Xj))/■\/Var(Xt) is 


expected to be distributed accordingly to a Gaussian distribution in the long-time limit. 
We plot on Fig. [T^the distribution V dehned by 


Vt{z) = Prob [Zt = z] , (142) 

and compare it with the normal distribution / a/^. At sufficiently long times, the 
distribution of the random variable Zt converges to the normal distribution of mean zero 
and unit variance. The rescaled position of the TP is then asymptotically Gaussian in 
the long-time limit. 
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Figure 10. Probability distribution function of the TP position Pt{x) = Prob[Xt = x], 
for different times. The results from numerical simulations (symbols) are compared to 
the expression (65), the functions w±i{u) being numerically computed as solutions of 
the system (124). The parameters are g = 0.15, pi = 0.98, p = 0.9. 
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Figure 11. Probability distribution of the rescaled variable Zt = [Xt — 
{Xt))l •\/Var(Xt) at different times, obtained with the following parameters: p = 0.1, 
g = 0.15, Pi = 0.98, T = r* = cr = 1. The black line is the normal distribution 
exp(—a;^/2)/v^. 
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7. Conclusion 

We studied the diffusion of a biased tracer particle (TP) in a hardcore lattice gas in 
contact with a reservoir of particles. From the general master equation of the problem, 
we gave a detailed derivation of the equation satished by the fluctuations of the position 
of the TP and presented in a previous publication. This equation involve the density 
prohles around the TP and cross-correlation functions, whose evolution equations are 
obtained in a closed form by resorting to mean-held type approximation. 

Going one step further, we extended this approximation to higher-order correlation 
functions in order to obtain the evolution equation verihed by the cumulant generating 
function of the TP position. This equation then yields the entire probability distribution 
of the TP on a lattice of arbitrary dimension. We also obtained the equation satished 
by the third-cumulant of the distribution, which gives information about its asymmetry. 

We then solved these equations in the particular case of a one-dimensional lattice. 
We recall the result from [21]: the dihusion coefficient of the TP is a nonmonotonic 
function of the density of bath particles. Counterintuitively, it reaches a maximum 
value for a nonzero value of the density. Thus, the presence of bath particles on the 
lattice may actually enhance its dihusion coefficient. Here, we showed that this ehect 
is related to an anomaly in the behavior of bath-tracer cross-correlation functions, that 
are nonmonotonic functions of the distance to the TP. Another surprising observation 
arises by analyzing the third cumulant of the TP position, which was shown to be 
nonmonotonic and to take negative values in a wide range of parameters. These 
analytical predictions were confronted with exact numerical samplings of the master 
equation, which indicate that the approximation we used is accurate in a wide range of 
parameters. 

We hnally solved the equation satished by the cumulant generating function, de¬ 
duced the probability distribution and showed that the position of the TP rescaled by 
its huctuations is Gaussian-distributed in the long-time limit. 

The equations presented in this paper are very general, and allow to compute the 
cumulants of the TP position and therefore its distribution under a mean-held-type 
approximation which was shown to be accurate in a broad range of parameters. These 
equations are valid for lattices of arbitrary dimension. We leave to future work the 
study of their solutions on higher-dimensional lattices. It could be interesting to see 
if the results obtained on the one-dimensional lattice - enhanced dihusion coefficient, 
nonmonotonic and negative third cumulant, convergence of the rescaled distribution to 
a Gaussian distribution - can be extended to lattices of higher dimension, or observed 
in experimental systems. 
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Appendix A. Evolution equations of {Xt^) 

In this appendix, we give an explicit derivation of Eq. (19), which governs the evolution 
of the second moment of Xt. We multiply the master equation (|^ by (X ■ eiY and 
average over all the bath conhgurations r] and all the positions of the TP X. We consider 
separately each term of the master equation: 

• the left-hand-side term of (|^ gives the contribution: 

C'l= 5^(X-ei)22dr*diP(X,77;f) 


x,v 
= 2dT*dt 


Y,iX ■ eifP(X,rr,t)] 

Kx,'n / 


= 2dT 


AX) 


dt 


(A.l) 

(A.2) 

(A.3) 


the hrst term of the right-hand-side yields 

C*! _ _ _ 

At=l XX-e^,X 


(A.4) 

(A.5) 


a 

= E [P{x,Ti"-^v)-P{x,mt)\ 

X,ri ^l=l XX-e^,X 

d 

= E(^'«itE E J2lX.ri''At)-P{X,r,;t)]. 

X /.=! r^X-e^,X V 

Recalling that a conhguration obtained from rj by exchanging the occupation 
variables of two neighboring sites r and r -|- e^, we obtain 

Y,P(Xpr;t) = Y.P(x,’i;t), 


(A.6) 


and we conclude that Ci = 0. 

the second term of the right-hand-side of Eq. (|^ gives the contribution: 

^(X ■ eif (1 - r]x) P(X - e^, p; t) 


Co 


■E 


Pv 


X,ri 


X,v 


(A.7) 


We consider for instance the term corresponding to /x = 1, and consider the hrst 
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sum over X and rj, in which we make the change of variable X X + ei: 

■ ei)2 (1 - T]x) P{X - ei, p-1) (A.8) 

x,v 

= Y.(X-ei+ a)^ (1 - ,,x+ej P{X,ri-, t) (A.9) 

X,v 

= 5^[(X ■ ei)' + + 2a(X ■ eO]' (1 - Vx+ejP{X, p; t). (A.IO) 

x,v 

Finally, we get 


■ ei)2 (1 - rix) P(X - e„i,i () - • e,)^ (1 - ',x+=,) P(X,rr, «)(A.ll) 

X,r] X,r] 

= 5^(1 - vx+er)PiX, p-1) + 2a ■ ei)(l - Vx+er)PiX, p-1). (A.12) 

X,r, X,r, 


Following the same procedure for the term p, = — 1 and noticing that the terms 
obtained for /i = ± 2 ,..., id in (A.7) cancel, we finally get 


9dT* 

C2 = - - keAt)) i 2 a((Xt) - geAt))] 

T 

+ P-i[(^^A - ke.AA) - 2 cx((Xt) - ge_At))]} , (A.13) 


where we define 

gr{t) = {XtT]Xt+r) ■ (A.14) 

• the third term of the right-hand-side of Eq. (|^ yields 

^3 =2dg 5^(X ■ eO' [(1 - Vr) P{X, 7 )- t) - PrP{X,p- f)] (A.15) 

X,r) r^X 

=2dg 5^(X ■ ei)2 [(1 - Vr) P(X, 77 ^ t) - VrP{X, 77 ; f)]. (A.16) 

X r^X r] 

Recalling that fj^ is the configuration obtained from 77 with the change Vr t— 1 — Vr^ 
we have the following equality 


(1 - Vr) P{X, 77^ t) = J2 VrP{X, 77 ; t), (A.17) 


which yields ^3 = 0 . 

• for the same reason, the fourth term will have a zero contribution after multiplying 
by (X ■ Ci)^ and averaging over X and 77 . 

Finally, bringing together the different contributions originating from the different terms 
of Eq. Q, we obtain 

- /^eiW) + 2a((Xt) -^eiW)] 

+p_i[a^(l - K_At)) - 2t^((^t) - de_i(t))]} , 


(A.18) 
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which is equivalent to Eq. (19): 


^ y {Pl [{^t) - C/ei {t)] -p_i[{Xt) - C/e_i {t)] } 


+ ^ {pi [1 - K, (t)] + p_i [1 - fce_i (t)] } 


(A.19) 


Appendix B. Evolution equations of Qr and Wr 

In this appendix, we start from the master equation (|^ in order to derive the 
evolution equations satisfied by gr{t) = {{Xt — {Xt))r]Xt+r) (Eq. (|24|) and Wr{u;t) = 
{D^^^rjxt+r) (Eq- (66)). The derivations of these two evolution equations are similar, 
and we will present a general method to derive the evolution equation of the following 
correlation function: 


frit) ^ {XiX,)px,+r) = J2j"(^)Vx+rPiX,rj;t), (B.l) 

X,v 

where is a generic function of the TP position. The equations satisfied by gr and Wr 
will be obtained by taking X{X) = X — (Xt) and X{X) = e™^ respectively. Paren¬ 
thetically, with this method, one retrieves the equation satisfied by krit) (Eq. ([^) by 
taking -T(A) = 1. 


We multiply the master equation (|^ by XiX)r]x+r and average over all the bath 
configurations rj and all the positions of the TP X. We consider separately each term 
of the master equation: 

• the left-hand-side term of (|^ gives the contribution: 

C-L = XiX)r]x+r2dT*dtPiX, p-1) (B.2) 

X,r, 

=2dT* PiX)px+rdtPiX, p- 1) (B.3) 

x,n 

= Vx+r [dt iPiX)PiX, p-t))- PiX, P- t)d,PiX)] (B.4) 

x,n 

In the two cases we will consider, we note that dtPiX) is independent of X and p, 
so that we obtain 


Cl =2dT*dt 


J2HX)vx+rPiX,p;t) 


lX,rj 

=2dT*dJrit) - 2dT*krit)dtPiX). 


- 2dT*dtPiX) Y, VX+rPiX, p- t) 

X,v 


(B.5) 

(B.6) 


where we defined fr in Eq. (B.l). We consider separately the two different possible 
expressions of P: 
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- if J^(X) = X-{Xt), then 


dtXiX) = 


d{Xt) 

dt ’ 


(B.7) 


and, using Eq. (§ , we obtain 

Cl = 2dT*dtgr{t) + 2dT*^ {pi [1 - K.it)] - p-i [l - K_,{t)] } (B.8) 

— if X{X) = e™^, then dtX{X) = 0, and we have the following expression of Cl: 

Cl = 2dT*dtWr{t). (B.9) 

the first term of the right-hand-side of the master equation becomes 


X fi=l r'^X-e^,X V 


P{X,p^'’^^;t)-P{X,p;t) . (B.IO) 


With an appropriate change of variable in the sum over rj, we obtain 

>?x+.-P(A, r,--'*; t) = Yi ’’’ 

V V 

so that Cl is written 


(B.ll) 


X ti=lr'j^x-e^,x n 

We consider the sum over r' 


V 


r ,fL 


X+r 


VX+r 


P{X,p;t). (B.12) 


E 


r’^X-e^,X 


X+r 


- Vx+T 


E 


V 


r ,/i 


X+r 


- VX+r 




X+r 


- VX+r 


(B.13) 


Recalling that a configuration obtained from rj by exchanging the occupation 
variables of two neighboring sites r and r -|- e^, we obtain the general relation 


7]:^ iir' ^ x,x - , 

{'n'' = { r]x-e^ if r' = X - , 

Hr' = x . 


(B.14) 


We then consider separately the two cases: 
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— if r = Cl, (z/ G {±1,..., ±(i}), using Eqs. (B.13) and (B.14), we obtain 


E 

r'i^X-e^,X 


r ,/! 


X+r 


- VX+r 


= '^fiVX+r + V_^? 7 x+r, 


(B.15) 


where the operator was dehned in the main text (Eq. ([^). We hnally 
obtain 

d 

Ci = E^wE {'^tiVx+r+ V-^r]x+r) P{X,r]]t) (B.16) 

X,rt At=l 

= 5^ {T{X)V,r,x+V 

V 

where the sum over /i runs over {±1 ,..., id}. 


(B.17) 

(B.18) 


— if r = Cl., using again Eqs. (B.13) and (B.14), we obtain 

d 


E E i’' 

^l=l r'^X-e^,X 

Then, Ci becomes 
6*1 = 


r ,/i 


X-\-e.u 


- VX+e^ 


V^?7x+e, - V-uVX+e.- (B.19) 


T(X) f^V^-V_„jr,x+, 


(B.20) 

(B.21) 


Finally, for any value of r, Eqs. (B.18) and (B.21) are recast under the equation 




(B.22) 


we then study the second term of the right-hand-side of the master equation (|^, 
which yields the contribution 


Co 


2dT* 


EE Pf,P{X)r]x+r 


X,r) fi 


X [(1 - Vx) P{X - e^, ri; t) - {l - 7]x+Y P{X, p; t)] 
2dT* 


(B.23) 


EE p,PiX)px+r (1 - Vx) P{X - e^, p-1) 


■ j-i 


(1 -hX+eJ PiX,p;t) l 
X,r, M J 


(B.24) 
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With the change of variable X ^ X + in the hrst sum, and recalling that 
X = X + Bi, we obtain 


^ 0 =- 


EE Pf,X{X + ae^ ■ ei)rix+r+e^ (l - Vx+e^) P{X, p-1) 


X,ri 


-^^Pt.P{X)px+r (1 -hX+eJ P{X,p;t) I . 
x,-n u J 

From this relation, we consider separately the different expressions of X: 
— in the situation where X(X) = X — (X^), we obtain 


(B.25) 


C'o=- 


- (Xt})px+r+e^ (l “ hX+eJ P{X, p' f) 


X,ri At 


'^^PviX - {Xt))px+r (1 - hX+eJ P{X, p- t) I 
X,ri At J 


'^'^PviX - {Xt))X^,px+r (l - hx+ej P{X,p;t) 


(B.26) 


X,t? A* 


+a EE ■ ei)px+r+e^ (1 - Px+e^) P{X, p]t)X 

X,r) 11 J 

As ■ Bi is equal to ±1 for ju = ±1 and 0 otherwise, we hnally obtain 

C 2 = —^pP(X, - (X,))(l - Vx.+eJ'^X.+r) 


(B.27) 


H - [Pl ((1 - PXt+ei)VXt+r+ei) “ P-1 ((1 “ PXt+e.i)PXt +r+e_i)] • 

(B.28) 

— in the situation where X{X) = e^"^^, we obtain 


C'o=- 


EE p,e“-^e‘“'’<'- -'i,x+.+.„ (1 - >)x+.„) P(X,,r,t) 


X,r) fi 


'^^PfiP'^^VX+r (1 -r/x+ej P{x,p;t) I 
X,r? ti J 


(B.29) 


(l - hx+ej - PX+r] P{X,p;t) 


X,r) fi 


(B.30) 


“ Vx+e^) V^Px+rP{X, p-, t) 


X,r] fi 

Orin-* ___ ___ 

- (l-hx+ej -l)px+r+e^P{X,p-,t). 


X,ri 11 


(B.31) 
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As ■ Cl is equal to ±1 for p, = ±1 and 0 otherwise, we finally obtain 

Orln-* __ 

+- - VXt+eJVX^+r+e,) • 

^ t=±l 


(B.32) 


the third term yields the contribution C^: 


C3=2dgJ2HX)vx.+r [{^-Vr')P{X,f^'-,t)-pr'P{X,v,t)\ (B.33) 

r'^X 


X r'^X 7J 


V 


X+r 


- VX+r 


P{X,p;t), 


(B.34) 


where we used again Eq. (B.ll). By the dehnition of fj'^' (conhguration obtained 
from Tj with the change rjrf 1 — Pr' ) , we write 


lix^ r', 

1 — Px if a; = r'. 

The expression of becomes 

C'a = 2dg ^ 7{X) ^ r]x+r{l - 2r]x+r)P{X, p; t). 

X V 

As px+r £ {0,1}, then {px+rY = Vx+r, and we hnally obtain 

Cs = -2dgfr{t) 


(B.35) 


(B.36) 


(B.37) 


the computation of the contribution (74 of the fourth term of the master equation 
is similar to that of C^, and we obtain 

Q = 2dfJ2HX)J2(^ - Vx+r){l -2vx+r)P{X,rj;t). (B.38) 


Using again {px+rY = Vx+r, we obtain 

(1 - hx+r)(l - 2T]x+r) = 1 “ VX+r, 

and, hnally. 


C, = 2df{{XiX,))-fA 
Note that when P{X) = X — {Xt), then {X{Xt)) = 0. 

Finally, writing from the master equation Q we obtain the relation 

C+ = Cl + C 2 + + C/i^ 

and we obtain the evolution equations satished by grit) and Wr(t): 


(B.39) 

(B.40) 


(B.41) 
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using (B.8), ( |B.22[ ), ( |B.28 ), ( B.37[ ) and ( B.40[ ), we get 
2dT*dtgr{t) = X] (V;, - gr{t) - 2d{f + g)gr{t) 


2dT* 


T 


2dT 


- VXt+e^)Vf,rix,+r) 


+- [pi +ei )vxt +r+ei) — P-1 a^-vx t+e_i )Vxt +r+e_i)] 


T 

2dT 

T 


-a 


{pi[l - K^{t)] -P-l[l - ke_,{t)]} kr{t). 


(B.42) 


which is equivalent to Eq. ( |24[ ), presented in the main text, 
using (B.9), ( B.22[ ), ( |B.32 ), (B.37) and (B.40), we get 


2dT*dtWr{u]t) = “ ^r,e^V_^^ Wr{u] 


]t) - 2d{f + g)wr (m; t) + 2df ) 


2dT* 


T 


2dT 


Pm - Vxt+e^)^^^VXt+r) 


Pe (e™"'" - 1) (e'“^‘(l - Px^+eJPXt+r+e,) 


(B.43) 


e=±l 


which is equivalent to Eq. ([6^, presented in the main text. 


Appendix C. Explicit expressions of gi and g^i in one dimension 

In this appendix, we solve the linear system of four equations satisfied by a, /?, and 


( 7 - 1 , obtained by substituting Eq. (95) into Eq. (93), Eq. (96) into Eq. (94) , and 


writing Eq. (95) for n = 1 and Eq. (96) for u = — 1. From Eq. (95) written with u = 1, 


and using the definition of W from Eq. (97), we obtain the following expression of a in 
terms of 'gi and 'g-p. 


a = 


Oq 


+ - + 


ai 


Aprx — A_\/r\ \ri Airi —A_i/ri 
where the coefficients Oq and a±i are defined by 

oo = A+- a 

T 


9i + 


a_i 


Airi - A_i/ri 


9-1, 


pi{l- p- K+ri){ri - 1) - p_i ( 1 - p - — ) (r^ - 1) 


K_ 


2D 


ai = 


a_i = 


T 

2D 

T 


PiK+{ri - 1 ) 


■p_iA+(ri - 1 ). 


(C.l) 

(C.2) 

(C.3) 

(C.4) 

(C.5) 


Similarly, from Eq. (|96|) written with n = — 1, and using the definition of W from Eq. 

isic 

hi 


(98), we obtain the following expression of ft in terms of and p_i: 
bo 




Air2-A_i/r2 Airs - A_i/r2 


9l+ [^2 


b-i 

Airs - A_i/rs 


9-1, 


(C.6) 
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where the coefficients bo and are dehned by 


2t* 


bo = - a 


T 


Plil - p - K+ri)(r2 - 1) -p_i ( 1 - p - — ) (r 2 ^ - 1) 


K_ 


hi = - 
b-i = 


2D 


T 

2D 


■piK_{r2 - 1 ) 


P-iK_{r^ - 1 ). 


In the boundary conditions (Eqs. (93) and (94)), replacing g 2 and p _2 by 


g2 = ari - 


2W 


-1 


ri 


AiVi - A_ir^ 

^13 2W' 1 

d-2 - 


r2^ 


Air2 - A^ir^ ^ r2^'' 


(C.7) 

(C. 8 ) 

(C.9) 

(C.IO) 

(C.ll) 

(C.12) 


and using the expression of a and {3 in terms of gi and p_i (Eqs. (C.l) and (C. 6 )), we 
obtain the following closed system of linear equations satished by and g-p. 


+ Bg^i — C 
Dgi + Eg_i = F 


(C.13) 


where we dehne 


2D 

A = - apiK+{ri - 1 )-- j—^ + Ain 

T Airi — A_ir^ 

A_i + 2{f + g) + —pi{p + K+rl) 

T 


2D 


B = - ap_iK+{r^ -1)- _ 

T Airi — A_ir;^ 


2 r* 

31 + —p-i{p + K+ri), 


2D 


C = - cr(pi(l - p- K+ri) - p_i(l - p- K_r2 )){p + K+ri) 

T 

-api(l - p- K+ri){p + K+rl) 


T 


Airf 


2D ( ( K \ K 

+- a ipi{l-p-K+ri)K+ri-p_i(l-p -^ . -d 

T \ V r2 J ri J Ain - A_in 


2D 


A ir 2 


-2 


D = - apiK_{r2-l) _ 

T Air2 - A_ir2 


31 + ‘^Pi{p + K.r^^), 


(C.14) 

(C.15) 


(C.16) 

(C.17) 


2D 


E = - ap-iK_{r2 

T 


-1 


1 ) ^-^^2 , A.r~^ 

Air 2 - A_ir 2 ^ ^ 

2D 


- {Ai + 2{f + g) + —P-i{p + K_r 2 ) 


(C.18) 
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and 


2t* 


F = -a(pi(l - p- K+ri) - p_i(l - p- K_r^ ))(p + K_r^ 

T 

+ ^ap_i(l - p- K_rf^){p + K_rf‘^) 


^-1^2 


-2 


2t* ( ( K_\ K_ 

- (T p-K+ri)K_r2-p-i(l-p -- . 

T \ \ r2 J r2 J Air2 - A_ir^ 


Solving the linear system given in Eq. (C.13), it is finally fonnd that: 

CE-BF 


9i = 


AE-BD 


and 9-1 = 


AF-CD 
AE-BD' 


(C.19) 


(C.20) 


K_, Ai and A_i have been determined in Section 5.1 Using their expressions, 


we compute gi and g-i, and finally compute K with the formula: 

2 

^ = ^ {pi{^ - P - K+ri) + p-i {l - p - K-rf^)) - ^ (pipi - p-ip-i). (C.21) 

Note that a and ft can be deduced straightforwardly from their relations with 'gi and 


g-i (Eqs. (C.l) and (C.6)) 


Appendix D. Explicit expressions of mi and m_i in one dimension 

In the main text we showed that, for n > 0, m„ had the following form 


mn = Trf + (a+n^ + b+n)rf, 


with 


a+ = 




2 Airi — A_ir^ ^ ’ 

K = ^-rT[C'i+ - a+(Airi + A_irf% 

Airi - A-ir^ 

and where the quantities Cf and Cf are defined by 

2t* _ _ 

Cf = - [pimiA+(ri - 1) +p_im_iA+(r]“^ - 1)] 


r 


2t* 


(D.l) 

(D.2) 

(D.3) 


- pia {2[(1 - A:i)(ari - aVi - a)- giK+fri - 1)] + cr(l - ki)K+{ri - 1)} 

r 

p_ia {-2[(1 - fc_i)(arf^ + a'rf'^ - a) - g-iK+{rf^ - 1 )] + (t(1 - k-i)K+{rf'^ - 1)} 


r 

2 r 


= — [pimiA’+(ri - 1 ) +p_im_iA + (ri ^ - 1 )] + V’+, 


(D,4) 


and 


Cf = - —cra'[pi(l - /ci)(l - ri) - p_i(l - ^)]. 

r 


(D,5) 


The expression of the coefficient a was given in Appendix C (Eq. (C.l)) and we 
defined a' = W/{Airi — A_iri“^). 
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For n < 0, a similar resolution leads to the expression 


with 


1 

a_ = -- 


b. = 


C-, 


1 = brf + (a_n^ + 6 _n)r 2 , 

(D. 6 ) 


(D.7) 

_i[Ci a_(Air 2 + A_ir 2 ^)], 

(D. 8 ) 


and where the quantities and C 2 are dehned by 

Cf = — [pimiK_{r2 - 1) + p-im-iK_{rf^ - 1)] 

Pia{2[(l - ki){(3r2- I3'r2 - /3) -giK_{r2 - 1)] + a(l - ki)K_{r2 - 1)} 
r 

-—p-io- {-2[(1 - k-i){/3rf^ + jd'rf^ - /3) - ^_iiF_(r 2 '^ - 1)] + cr(l - k-i)K_{rf^ - 1)} 


r 

2 r 


= — [pimiiF_(r2 - 1) +p_im_iiF_(r2 ^ - 1)] 


(D.9) 

(D.IO) 


and 


^ 2 - = 


AD 


T 


a^'[pi(l - fci)(l - rs) - p_i(l - - r 2 % 


(D.ll) 


The expression of the coefficient ft was given in Appendix C (Eq. (C. 6 )) and we 
dehned ft' = Wj{Apr 2 — A_ir 2 “^). 

As explained in the main text, we hnally a linear system of four equations with 


unknowns mi, m_i, F and b (Eq. (Ill)), involving the following quantities: 

2t* 


Mi3 — 


ri 


T Airi — A_ir^ 


^\piK+{ri - 1) - 1] 


Ml 4 — 


2t* 


ri 


M 23 — 


T Airi — A_ir^ 
2t* 1 


^PiK+{r^ - 1 ) 


r r2 Air2 - A_ir2 


3YPiiF_(r2 - 1) 


Moa — 


2t* 1 


r r 2 Air 2 - A_ir 2 


^[P-iK-{r2 - 1 ) + 1 ] 


(D.12) 

(D.13) 

(D.14) 

(D.15) 


2 Ai 


2t* \ 2t* 

M 33 = — ( A_i + 2{f + g) - pik 2 ) H- - 

T J T AiVi — A-iV]^ 


3 YPiiF+(ri - 1) (D.16) 


2r* , 2t* 2Ai 

Af34 = - p-iki H---- _ 

T T Airi — A_ir;^ 

2 A i 


^P-iK+{r^ - 1 ) 


,. 2 r* , 2t* 

Mas = - P-iki -, , 

T r 4ir2 — A-ir., 


—PiA'_(r 2 - 1 ) 


(D,17) 

(D.18) 















Distribution of the position of a driven tracer in a hardcore lattice gas 


48 


A 444 — 


-(A, + 2{f + g) + _ 

T J 'T Mr2 - ^-1^2 


ziiP-iK-ir^ 


-1 


>"2 = - 
^3 = —pi — 2^1 


>4 = —(^-1 — 2 ^ 4 . 


2a+y4_i — fjj^ri 
Airi - A_irf^ 

2a_Ai — 

Airs - A_irf^ 
a+{Airi - 3A_irj~^) + 

Airi - A_irf^ 
a_(3Air2 - A_ir^^) + 'ijj_ 


Air2 - A_ir2 


-1 


1) (D.19) 
(D.20) 
(D.21) 
(D.22) 
(D.23) 


Appendix E. Equations satisfied by w±i{u): additional definitions 


In this appendix, we give the explicit expressions of the coefficients Bj, Cj, Dj, Ej 
involved in Eq. (124). 


i?i = Ai + —pi(e‘"'^-l)(l-fci) 

T 

-fs 3fs 

B 2 = Ai+A_i+ 2 (/+ 5 () + —pi(e'“'^-l)[l-w;i(M)] + —p_i(e"'“^-l)[l-w;_i(M)] 

T T 

B, = A_^ + - 1)(1 - fc_i) 

T 

S4 = ^ {pi[ki - wi{u)]{rie'^°' - 1) +p_i[fc_i - M;_i(M)](r]“^e“'“'" - 1)} 

^5 = 2/ + {pi[ki - tri(M)](e'“'" - 1) +p_i[fc_i - w;_i(M)](e"'“'" - 1)} 

<^4 = ^ {pi[ki - Wi(M)](r2e™'" - 1) +p_i[/c_i - w;_i(M)](r2'^e“'“'^ - 1)} 

Di = Ai + —pi(e™"-l)(l-A;i) 

T 

D 2 = A_i + 2(/ + g) + —pie‘""fc 2 + — [pi{E^^ - 1) + p_i(e-'“'^ - 1)] 

TT 

2t* 

D 3 = - P-iki 

T 

E >4 = 2 / + ^ [pie™'"/cifc 2 - P-ikik_i\ 

Ds = —Pi(e™" - 1) 
r 

De = —p_i(e-“" - 1) 

T 

r. 2r* , 

El = -pifc_i 

T 


(E.l) 

(E.2) 

(E.3) 

(E.4) 

(E.5) 

(E.6) 

(E.7) 

(E.8) 

(E.9) 

(E.IO) 

(E.ll) 

(E.12) 

(E.13) 










Distribution of the position of a driven tracer in a hardcore lattice gas 


49 


E 2 = A^ + 2{f + g) + —p_,e-^--k .2 + — - 1) + - 1)] (E.14) 

T T 

Es = A_i + —- 1)(1 - fc_i) (E.15) 

r 

£’4 = 2/+ — \p_ie~''^'"k-ik-2 -pikik_i] (E.16) 

r 

Appendix F. Algorithm and numerical methods 

To sample exactly the master equation Q, we generate a sequence of random numbers 
{T,x,fi) with the joint probability density function p(r, x, p) where p(r, a;,/r)dr is the 
probability at time t that the next event occurs in the inhnitesimal time interval 
[t + T,t + T + dr], at site x, and is of type fv (i.e. a diffusion event, an absorption 
event or a desorption event). We write 

p ( r , x ,/ i ) = pi{t)p 2 {x\t)p 3 {h\,x,t) ( F . l ) 

where 

• pi(r)dr is the probability at time t that the next event occurs in the time interval 
[t + T,t + T + dr]. 

• P 2 {x\T)dT is the probability that the next event occurs at site x, knowing that it 
occurs during the time interval [t + T,t + t + dr]. 

• p 3 (p|x,r)dr is the probability that the next event is of type /i knowing that it 
occurs during the time interval [t + T,t + t + dr] and at site x. 

Writing Cx,i_i the transition rate of event /i at site x, we dehne 

V 

R=Y,rx. ( F . 3 ) 

X 

Vx is then the total rate of the events at site x, and R is the total rate of all the events 
on all lattice sites. The quantities r, x and pi are then respectively drawn from the 
following distributions: 


PIW 

= Re-^\ 

(F.4) 

P2{,x\t) 

Tx 

~ 

(F.5) 

P3(hk,r) 

ir> 

(F.6) 


The algorithm is as follows. We build a lattice of length 2L + 1 (the spacing of the 
lattice a is taken equal to 1). The boundary conditions are periodic, and L is chosen to 
be large enough so that we can consider the lattice as inhnite (in the results presented 
below, L > 250). The initial condition is the following: the TP is initially at the origin 
and at each site different from the origin, a particle is set with probability p. We chose 
a hnal simulation time tmax- At each step of the simulation, and as long as t < tmax, the 
algorithm follows these steps: 
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(i) Set R = 0. 

ii) For each x G [0, 2L + 1], compute r^: 

• if the site x is occupied by a bath particle, three events are possible : a jump 
to the left, a jump to the right, or a desorption event (respectively labeled by 
1, 2 and 3). The associate rates are 
1 


111 


vni 


Cx,i = -(1 - g)(l -rjx-i) 
Cx,2 = ^(l-g)(l-T]x+i) 


(F.7) 

(F.8) 

(F.9) 


^x,3 9 

SO that = ^(1 - g){l- Px-i) + ^(1 - - Px+i) + 9- 

• if the site x is occupied by the TP, two events are possible : a jump to the left, or 
a jump to the right. The local rate is then r,j, = (1—pi)(l—pa,_i)+pi(l —pj.+i). 

• if the site is empty, the only possible event is an absorption event, and r^ = /. 

Compute the total rate R = 

(iv) Draw r from the distribution Pi(t) = Re~^'^. 

(v) Draw x from the distribution P 2 {x\t) = r^/R. 

(vi) Draw the event pi from the distribution p 3 (p|x,r) = Cx,^/rx- 

vii) Update the lattice occupation after the realization of the event (note that such an 
algorithm is rejection-free). 

Increase the time t t + r. 

With a large number of 


We hnally keep track of the TP position Xt with time 
realizations, we sample the p.d.f. of Xt. 


Appendix G. Range of parameters for which K > 1/2 


In section 6.3.3 we determined the critical value of the desorption rate gc allowing the 
emergence of a nontrivial maximum value for the function K{p) for a given value of 
the bias 6. As a complementary approach, for a hxed value of the bias, we study the 
domain of the plane {9,p) in which the diffusion coefficient is greater than 1/2 , which 
is its value when there is no bath particle. 

In the case where the curve K{p) displays an extremum value, the value of the 
diffusion coefficient is greater than 1/2 in the range [0,p], where p is a function of pi 
and g only. From the analytical curves K{p), we can then deduce the values of p as a 


function of g for different values of pi, which are plotted on Fig. G1 


As it was shown on Fig. for a given value of pi and p, there exists a value of p 
below which the difference p_i — p _2 is negative, i.e. below which there is a minimum 
for 'gn in the range n < 0. This critical value of p will be denoted by /f. In order to 
compare the domain of parameters giving respectively K > 1/2 and p_i — p _2 > 0, we 


also represent on Fig. G1 the curves of p as a function of g for different values of pi. 





Distribution of the position of a driven tracer in a hardcore lattice gas 


51 
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9 

Figure Gl. Critical value of the density p below which the diffusion coefficient is 
greater than 1/2, as a function of g and for different values of the bias pi (curves in 
blue). Critical value of the density below which — g _2 > 0, as a function of g 
and for different values of the bias pi (curves in green). 


The domains in the plane {p-,g) for which K > 1/2 and for which g^i — g -2 < 0 
then seem to be correlated. 
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